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Abstract 

We developed a prototype time projection chamber using gas electron multipliers 
(GEM-TPC) for high energy heavy ion collision experiments. To investigate its 
performance, we conducted a beam test with 3 kinds of gases (Ar(90%)-CH4(10%), 
Ar(70%)-C2Hg(30%) and CF4). Detection efficiency of 99%, and spatial resolution of 
79 /im in the pad-row direction and 313 /im in the drift direction were achieved. The 
test results show that the GEM-TPC meets the requirements for high energy heavy 
ion collision experiments. The configuration and performance of the GEM-TPC are 
described. 
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1 Introduction 



High particle multiplicity is an important feature to be considered in de- 
tector designs for high energy heavy ion collision experiments. In ^/snn = 
200 GeV Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory, the average charged particle multiplic- 
ity (dN ch /drj\ v=0 ) is 170 Therefore, average charged particle density is 
0.03 cm -2 at a distance of 30 cm from the vertex. Additionally, experiments 
at RHIC are performed at a high event rate of about 10 kHz and the charged 
particle rate is 300 cps/cm 2 at a distance of 30 cm from the vertex. This harsh 
environment demands highly efficient central tracking detectors. 

A wide variety of observables are measured in high energy heavy ion collision 
experiments: Such as charged particle multiplicities, yield ratios and spectra 
of identified hadrons, elliptic flow, suppression of high p T particle production 
and heavy flavor production 0. Therefore, several particular features are re- 
quired for the detectors used. A relatively wide transverse momentum range 
(px ~ 0.2-20 GeV/c) is required to be covered to take in the broad interests 
of high energy heavy ion collisions and the magnetic field should be kept low 
(~1 T). However, good momentum resolution of 5pt/pt ~ 3x 10~ 3 (GeV/c) -1 
is required for future measurements, such as the T state measurement at 
RHIC-II [3J. To achieve such good momentum resolution with a magnetic field 
of 1 T, a 1-m radius solenoidal tracker with spatial resolution of ~200 fim is 
needed. Double track resolution is also required to be better than 1 cm to cope 
with the high particle multiplicity. Such high flux operational and multi hit 
capabilities have recently also been required in particle physics experiments 0]. 

One of the candidates as a detector that satisfies the above requirements is 
a combination of a time projection chamber (TPC), a sophisticated detector 
for particle tracking and particle identification, and micro patter gas detectors 
(MPGDs) jEf such as gas electron multipliers (GEMs) jfjj. Although existing 
wire chambers need gating wires to collect positive ions and have limitations 
in their double track resolution by the wire spacing, the novel structure of a 
GEM has the following advantages in its application to TPC readout: 

• Two dimensional symmetry of the GEM structure provides a large flexibility 
in the readout geometry. 

• Large areas can be covered with a low amount of material since support 
mechanics for a GEM are simple. 

• The intense electric field region will be limited inside the holes and the 
ExB effect is expected to be reduced. (In a strong magnetic field, this 
effect results in a broadening of the electron cloud and a worsening of the 
resolution.) 

• The signal of the readout pad is induced by amplified electrons, and is 
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spatially narrow and fast. 

• The generated ions do not affect the signal and the ion tail of the signal will 
be suppressed. 

• The positive ion feedback into the drift region can be suppressed by a factor 
of about 10 due to the electric field around the GEM holes and so gating 
wires might be unnecessary. 

Therefore, a TPC using a GEM for signal amplification (GEM-TPC) may 
achieve high rate capability as well as excellent double track and spatial res- 
olution. A GEM-TPC is a strong candidate to be a central tracking detector 
in high energy heavy ion collision experiments 0, H, 0, Ell II • 



2 GEM-TPC Prototype 



2.1 Mechanical Structure 



A GEM-TPC prototype, consisting of an end cap chamber, a gas vessel and 
a field cage, was developed [12| • Photographs of the GEM-TPC prototype are 
shown in Fig. 1. Figure 2 shows a schematic view of the GEM-TPC. 

The end cap chamber mounts a triple GEM (the effective area is 10 x 10 cm 2 ) 
on readout pads. The triple GEM was constructed from GEM foils made at 
CERN As shown in Fig. 2, the gap between neighboring GEMs was 
2 mm and the gap between the bottom GEM and the pads was 1.5 mm. High 
voltages are applied to the triple GEM through connectors penetrating the 
end cap chamber with a resistor chain. When the voltage across the GEMs 
(Vgem) is 360 V, the electric fields in the transfer region and the induction 
region are ^=360 V/2 mm=1.8 kV/cm and -Ej=360 V/1.5 mm=2.4 kV/cm, 
respectively. The operated drift electric field is shown in Table 1. Two kinds 
of readout pads with different shapes, rectangle and chevron (zigzag), were 
used to study the dependence of the spatial resolution on shape (see Fig. 3). 
Since chevron pads may increase the number of hits on pads by charge sharing, 
chevron pads are expected to have better spatial resolution than rectangular 



ones 



14j . Both kinds of pads, both made of gold-plated copper, have the same 
area of 1.09 x 12.0 mm 2 and the same pitch of 1.27 mm. Relatively narrow 
width pads are required for charge sharing for small diffusion gases such as 

cf 4 O US 

The field cage is a cuboid with dimensions of 36 (drift direction) x 17 x 17 cm 3 
and creates a uniform electric drift field. The electric field uniformity is 
Jq c ™(E±/ E)dz < 1 mm in the center area of 10 x 10 cm 2 which corresponds 
to the GEM effective area. The field cage consists of 115 gold-plated copper 
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Fig. 1. The gas vessel of the GEM-TPC (top left), the field cage (top right) and the 
end cap chamber (bottom). 



strips connected with 1-MQ resistors in series on FR4 boards. At the end of 
the resistor chain, additional resistors are placed to match the voltage at the 
bottom of the field cage with the surface voltage of the top GEM. 

The gas vessel is made of aluminum plates and has a volume of 60 x 29 x 29 cm 3 . 
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Fig. 2. A schematic view of the GEM-TP C. 
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Fig. 3. Readout pad layout. There are 3 rows x 10 columns for both rectangular 
and chevron pads. The outer 2 columns are not read out. Each pad has an area of 
1.09x12.0 mm 2 . 



2.2 Front End Electronics 



A charge sensitive preamplifier, consisting of two kinds of operational ampli- 
fiers (AD8058 and AD8132, Analog Devices, Inc.), was used for the GEM- 
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Table 1 

Properties of the 3 kinds of gases. Drift velocities and diffusion coefficients were 
calculated by Magboltz. 

Gas Ar(90%)-CH 4 (10%) Ar(70%)-C 2 H 6 (30%) CF 4 

Operated drift field (V/cm) 130 390 570 

Drift velocity (cm//is) 5.48 5.01 8.90 

Transverse diffusion (fim/^cm) 570 306 104 

Longitudinal diffusion (/xm/^/cm) 378 195 82 

Mean energy for ion-electron 26 26 54 

pair production (eV) 

TPC. Its time constant is r = 1 fis and its effective gain G = 3.2 V/pC. The 
gain of the preamplifier was determined from the expected signal amplitude 
and the dynamic range (0 ~ -1 V) of a flash ADC (FADC) module (RPV- 
160, REPIC Co., Ltd.). The resolution and the sampling rate of the FADC are 
8 bits and 100 MHz, respectively. Signals from 24 pads (3 rows x 8 columns) 
are transmitted from the preamplifiers to the FADCs through 8-m shielded 
twisted cables. 



2.3 Gas 



Three kinds of gases with different properties, a mixture of argon(Ar)(90%)- 
methane(CH 4 )(10%) (commonly called P10), a mixture of argon(Ar)(70%)- 
ethane(C 2 H 6 )(30%) and pure tetrafluoromethane(CF 4 ) (99.999%) gas, were 
used to study the performance of the GEM-TPC. Properties of these gases 
are shown in Table 1. Drift velocities and diffusion coefficients were calculated 
by Magboltz Q. 

Ar(90%)-CH 4 (10%) is widely used in TPCs and GEMs jLH, El Q. A fast 
drift velocity that peaks at a low electric field of 130 V/cm is the primary 
attribute of Ar-CH 4 , and such a low electric field is an advantage for de- 
tector operation. However, its large diffusion coefficients are a disadvantage. 
Ar(70%)-C2H 6 (30%) is a common chamber gas, but there are few results for 
it with a GEM. Although its diffusion coefficients are smaller than those of 
Ar-CH 4 , the drift velocity peaks at a relatively high electric field of 550 V/cm. 
CF 4 gas is studied as a TPC and GEM gas because of its very small diffusion 
coefficients and very fast drift velocity (l^, EH . These properties of CF 4 will 
be advantages in high particle multiplicity environments of high energy heavy 
ion collision experiments. However, a high electric field is needed to achieve a 
fast drift velocity with CF 4 . 
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The gas flow rate was set at ~200 ml/min using a mass flow controller (SEC- 
E40, ESTEC Co., Ltd.). Gas pressure was set at the atmospheric pressure 
using a bubbler filled with silicone oil. 

Effective gas gains were measured using an 55 Fe X-ray (5.9 keV) source. An 
X ray creates primary electrons by a photoelectric effect with a gas molecule. 
The obtained 55 Fe X-ray charge spectrum with CF 4 is shown in Fig. 4. The 
spectrum around the 5.9-keV peak is fitted by a Gaussian function and the 
gain G is calculated as 



G _(Q) W 

where (Q) stands for the mean of the Gaussian function, e is the electron 
charge, and W is the mean energy for ion-electron pair production and is 
given in Table 1. The gains for the 3 kinds of gases are shown as a function of 
Vqem in Fig. 5 and a gain of 10 4 was achieved with each gas. Although CF 4 
needs a higher voltage than Ar-CH 4 and Ar-C2H 6 to achieve the same gain, 
the slopes of the gain curves are similar for the 3 gases. From the obtained 
Gaussian sigma, erg, the relative energy resolution (a E /E) with Ar-CH 4 , Ar- 
C 2 H 6 and CF 4 is determined to be 11%, 10% and 13%, respectively. The 
obtained energy resolution is comparable with the result by another research 
group (16% for CF 4 and 9% for Ar(70%)-CO 2 (30%)) 0. 



3 Performance Test 



A beam test was performed at the ir2 secondary beam line of the 12-GeV 
Proton Synchrotron at KEK (KEK-PS) to evaluate the basic performance of 
the GEM-TPC. 

The characteristics of the GEM-TPC evaluated in the performance test were 
detection efficiency, spatial resolution in the pad-row direction and the drift di- 
rection and particle identification capability by dE/dx measurement. The de- 
pendence of these characteristics on the 3 kinds of gases (Ar(90%)-CH 4 (10%), 
Ar(70%)-C 2 H 6 (30%) and CF 4 ), the GEM gain (6x 10 2 -2x 10 4 ), the drift length 
(20-290 mm), the readout pad shape (rectangle and chevron), the beam mo- 
mentum (0.5-3.0 GeV/c) and the beam rate was also evaluated. Except for 
the GEM gain dependence measurement, Vqem was fixed to 332 V, 341 V and 
498 V for Ar-CH 4 , Ar-C 2 H 6 and CF 4 , respectively. 

Figure 6 shows a schematic view of the detector setup for the performance 
test. Three plastic scintillation counters (SI, S2 and S3) were used for event 
triggering and time of flight measurements for particle identification. Two 
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Fig. 4. Measured 55 Fe X-ray (5.9 keV) charge spectrum with CF4. The energy 
resolution is <je/E = 13%. 

gas Cherenkov counters (GCC1 and GCC2) filled with 2.5-atm CO2 gas were 
used for electron identification. Two silicon strips detectors (SSD1 and SSD2), 
each of which consists of two single sided strips (strip pitch of 80 jum) crossing 
at right angles, were used for particle tracking. Two hadron blind detectors 
(HBD1 and HBD2) were tested at the same time and the results are shown 
in [20|. The GEM-TPC was operated without a magnetic field. 

Figure 7 shows a typical signal of the GEM-TPC operated with Ar-C2H 6 . The 
signal was recorded for 6.4 fxs (=640 samples) in one event. One channel of 
the ADC corresponds to ~ 4 mV. To extract the pulse height and the arrival 
time of the signal, the following function is fitted to the FADC spectrum, 



ADCn-n + n-^-f-rt™ (2) 

1 + exp(-(t -P2J/P4) 

where t is sampling time. The fitting parameters in the above function can 
be recognized as follows: po is the pedestal, p\ is the pulse height, P2 is the 
arrival time, p% is the time constant of the electronics and p^ is the rise time. 
The obtained pulse height and arrival time are used for determinations of the 
hit position. The hit position in each pad row is determined by an amplitude 
weighted mean of pad positions and arrival times: 
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Fig. 6. A schematic view of the detector setup in the performance test. 



Xi 



in the pad-row direction and 



(3) 



_ 2^=0 Pl,i,j ■ P2,i,j ■ V drift ( .s 

Zi ~ V 7 Th ■ ■ { ' 

2^j=oPl,i,J 

in the drift direction, where Pk,i,j is the k-th parameter of the i-th row and the 
j-th column pad (0 < i < 3 and < j < 8), D = 1.27 mm is the pad spacing 
(see Fig. 3) and v drift is the drift velocity. 
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Fig. 7. A typical GEM-TPC signal crated by a 1-GeV/c electron beam. The drift 
gas is Ar(70%)-C2He(30%), the drift length is 85 mm and the pad shape is chevron. 
A voltage across the GEMs is Vqem = 341 V. 



4 Results 



4-1 Detection Efficiency 



Single-pad-row detection efficiency was measured as a function of Vqem with 
the 3 kinds of gases. Measurements were done with 1-GeV/c ir~ beams with 
a drift length of 20 mm with Ar-CH 4 and 85 mm with both Ar-C2H 6 and 
CF4. Tracks having hits in the 1st and 3rd pad rows were selected for the 
efficiency evaluation, and the fraction of the hits in the 2nd pad row was 
used as the detection efficiency. Results are shown in Fig. 8. The detection 
efficiency reaches a plateau at a gain of ~ 4 x 10 3 with Ar-CH 4 and Ar- 
C2H6. The small diffusion of CF4 makes the efficiency reach a plateau at a 
smaller gain of ~ 2 x 10 3 . The efficiency plateaus are 99.3%, 99.6% and 99.8% 
with Ar-CH4, Ar-C2Hg and CF4, respectively. These results are very similar 
to results from another research group findings of Ar(70%)-CO 2 (30%) and 
Ar(93%)-CH 4 (5%)-C0 2 (2%) 0. 

The gain scales for Ar-CH 4 and Ar-C2H 6 in Fig. 8 are determined by the 
obtained gam-Vc em relation in Fig. 5 and the one for CF 4 is determined by 
measured pulse heights using the relation between the gains and measured 
pulse heights with Ar-CH 4 and Ar-C2H 6 . The gain with CF 4 in Fig. 8 is 
larger than the one in Fig. 5. This difference may be due to differences of 
the atmospheric pressure p and the temperature T because the GEM gain 
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strongly depends on p/T 19j|. Unfortunately, the atmospheric pressure and 
the temperature were not monitored during the performance test. 
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Fig. 8. Detection efficiency of the GEM-TPC as a function of a GEM gain. 



4-2 Transverse Diffusion Coefficient 



Transverse diffusion coefficients were measured with 1-GeV/c tt~ beams. The 
coefficients are evaluated using the spatial distribution of secondary electrons 
in the pad-row direction. The secondary electron distribution is fitted by a 
Gaussian distribution. The obtained sigma of Gaussian, s x , can be expressed 
as 



si(L) = si + C 2 DT ■ L } (5) 

where L is the drift length, Cbt is the transverse diffusion coefficient and 
s x q is the intrinsic width of the induced charge distribution determined by 
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Table 2 

The measured width of the induced charge distribution and transverse coefficients. 
Only statistical errors are shown. 



Gas 


Pad shape 


SxO 


(jim) 


C DT (H 


\Jcm) 


Ar(90%)-CH 4 (10%) 


chevron 


385 


± 101 


508 ± 


7 


Ar(90%)-CH 4 (10%) 


rectangle 


387 


± 101 


505 ± 


7 


Ar(70%)-C 2 H 6 (30%) 


rectangle 


402 


± 43 


317 ± 


4 


CF 4 


chevron 


383 


± 37 


107 ± 


6 



the readout system configuration. Figure 9 shows si as a function of the drift 
length. The measured s x0 and Cdt are shown in Table 2. The measured values 
of Cdt of Ar-C2H 6 and CF 4 agree well with the calculated values of Cdt shown 
in Table 1; however, for Ar-CH 4 the measured one is slightly smaller than the 
calculated one. 
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Fig. 9. The squared width of the diffused secondary electrons s% for the 3 kinds of 
gases. 



4-3 Spatial Resolution 



Single-pad-row spatial resolution in the pad-row and drift directions was eval- 
uated for a drift length range of 20-290 mm with 1-GeV/c 7r~ beams. The 
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Table 3 

The number of effective secondary electrons, N e ff, and the number of secondary 
electrons, N, for 1-GeV/c n~ beams and the track length of 12 mm. 



Gas Pad shape N eff N N eff /N 



Ar(90%)-CH 4 (10%) 


chevron 


31± 1 


119 


0.26±0.01 


Ar(90%)-CH 4 (10%) 


rectangle 


30± 1 


119 


0.25±0.01 


Ar(70%)-C 2 H 6 (30%) 


rectangle 


31± 1 


131 


0.23±0.01 


CF 4 


chevron 


86±23 


147 


0.58±0.15 



single-pad-row spatial resolution in both directions was evaluated by the resid- 
ual between the hit position of the 2nd pad row and the interpolated hit posi- 
tion from the 1st and 3rd pad rows. The measured spatial resolution is shown 
in Fig. 10. The effect of diffusion on the spatial resolution in both directions is 
clearly seen. The best resolution is 79 /im in the pad-row direction and 313 fim 
in the drift direction, obtained with Ar-C2H 6 gas and rectangular pads at 13- 
mm drift. The spatial resolution of the chevron pads is almost the same as 
that of the rectangular ones. If the charge distribution is Gaussian, the spatial 
resolution of the chevron pads should be better than that of the rectangular 
] . A possible reason for the result obtained is that the finite sizes of the 
GEM holes distort the charge distribution from a Gaussian distribution and 
so the non-Gaussian tails worsen the spatial resolution of the chevron pads. 

The dependence of the spatial resolution in the pad-row direction on the drift 
length of L can naively be understood as 



<(L) = v'o + CbT-L/Neff, (6) 

where a x0 is the extrapolated resolution at a zero drift length and N e jf is the 
effective number of secondary electrons 11 1. N e ff is mainly determined by 



the number of secondary electrons, N, and pad geometry [21]. The calculated 
number of secondary electrons, N, and the measured number of effective sec- 
ondary electrons, N e ff, for 1-GeV/c tt~ beams and a track length of 12 mm 
are shown in Table 3. A small fraction of the secondary electrons effectively 
contribute to spatial resolution. 



4-4 Beam Rate Dependence 



One of the advantages of the GEM-TPC is its ion feedback suppression. The 
effect of ion feedback on GEM-TPC performance was studied by measuring 
the beam rate dependence of the detection efficiency and spatial resolution. 
The beam rate was determined by the beam slit width and the rate was 
monitored with the S2 scintillator (2.5x2.5 cm 2 ). The beams of e + , 7r + andp at 
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Fig. 10. Spatial resolution in the pad-row direction (top) and the drift direction 
(bottom). 

a momentum of 2 GeV/c and Ar-CH 4 gas were used, while the drift length was 
85 mm. The results are shown in Fig. 11. The results of the detection efficiency 
and the spatial resolution in the previous subsections were obtained with a 
beam rate of < 500 cps/cm 2 . At the maximum beam rate of 4800 cps/cm 2 , 
the detection efficiency and the spatial resolution were worsened by factors 
of 2.5±0.5% and 11±3%, respectively. The maximum total beam rate in the 
active GEM-TPC region was in the order of 10 5 cps. 

This result is worse than the results from the other research group with beams 
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Fig. 11. Dependence of the detection efficiency (top) and the spatial resolution on 
the beam rate. 

of 1.5 x 10 5 cps/cm 2 0| . Because of the limited readout area of our GEM-TP C 
(1x3.6 cm 2 ), it was not possible to fully distinguish double tracks, which 
worsened the detection efficiency and the spatial resolution. 

The beam rate exceeds the typical rates in y / s/v TT =200 GeV Au+Au collisions 
at RHIC and v /sjvat=5.5 TeV Pb+Pb collisions at LHC, which are 300 cps/cm 2 
and 1400 cps/cm 2 , respectively, at a distance of 30 cm from the vertex. Since 
4800 cps/cm 2 is much larger than these numbers, the effect of the ion feedback 
on the GEM-TPC performance can be regarded as negligible for our purpose. 



15 



4-5 Particle Identification Using dE/dx 



Energy losses, dE/dx, were measured for positrons, muons, pions, protons 
and deuterons in a beam momentum range of 0.5-3.0 GeV/c to evaluate the 
particle identification capability. In this measurement, the drift length was 
85 mm and Ar-CH4 was used. The summation of pulse heights for 3 pad 
rows is regarded as the energy loss. Since there was a time variation in the 
GEM gain (30% peak to peak, for 5 hours), the measured value of the mean 
pulse height, (PH) l meas (i = e + , 7f + ,p, d), was corrected with pions at each 
momentum to eliminate the effect of this time variation as follows. 



(PH) n+ 

(PHY = (PHY ■ - '-2*. (7) 

\ I corr \ I meas i p TT\ 7r + ' \ / 

\ / meas 

where (PH) l corr is the corrected mean pulse height and (PH) l exp is the expected 
mean pulse height. Figure 12 shows the corrected mean pulse heights and 
curves of the expected mean pulse heights. As previously mentioned, the gain 
variation may be due to the change of p/T. Another candidate as a reason for 
the gain variation is the charging- up of the insulator of a GEM foil |2^ |. 




q LJ i i i i i i i i_l i i i i i i i l_L 

10~ 1 1 10 

Momentum [GeV/c] 

Fig. 12. Corrected mean pulse heights for 5 kinds of particle species. The curves are 
the expected mean pulse heights. 

To estimate the particle identification capability of a large GEM-TPC, a 
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Monte-Carlo simulation was performed using measured pulse height spectra 
for 1-GeV/c pions and protons. To improve the energy resolution, a truncated 
mean method, where 2/7 of the pad rows having the largest signals are re- 
moved, was used. Energy resolution of pions will be 9% and the pion rejection 
factor with 99% proton efficiency is expected to be 200 with a 50-cm track 
length. This energy resolution is comparable with that of the STAR TPC (8%) 
with a track length of more than 67 cm (l7| . 



5 Conclusion 

A GEM-TPC prototype was constructed to develop a tracking detector for use 
in high event rate and high particle multiplicity environments in high energy 
heavy ion collisions. 

To evaluate the performance of the GEM-TPC, a beam test was performed 
at KEK. Detection efficiency of >99.3% was achieved with 3 kinds of gases, 
Ar(90%)-CH 4 (10%), Ar(70%)-C 2 H 6 (30%) and CF 4 . Spatial resolution of 79 /an 
in the pad row direction and 313 /tm in the drift direction was achieved 
with Ar-C2Hg and rectangular pads for 13-mm drift. The GEM-TPC showed 
high detection efficiency and good spatial resolution with a particle rate of 
4800 cps/cm 2 , which exceeds the particle rate of RHIC and the LHC. Energy 
loss measurements showed a good particle identification capability. 

These results indicate that the GEM-TPC meets the requirements for central 
tracking detectors for use in the next generation of high energy heavy ion 
collision experiments. 
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